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^ • Abstract 

The exact form of the form-factor associated with a meson-photon-photon 
vertex is not known. Different suggestions exist, based on constraints from 
' QCD and on recent experiments. Four different form-factors are studied in 

^ this article. 

We calculate decay rates for the decays 7r° — > 77, 77 7 e+e~, 
7] — fi~^fi~ e~^e~ and 77 e"*"e~ e"'"e~ and cross sections for the processes 
e+e~ PS e+e~, where PS = ir^, rj or rj'. The results depend on the 
choice of form-factor and we examine if the differences are large enough to 
distinguish in upcoming experiments. 

The dominant part of the light-by-light contribution to the muon anomalous 
magnetic moment is also calculated. Here the uncertainty around the choice 
of form-factor implies that it is possible to estimate a smallest error in the 
theoretical value of this contribution. 



^Master of Science Thesis by Fredrik Persson with thesis advisor: Johan Bijnens. A more paper hke 
version with more references will follow. 



1 Introduction 

Understanding the structure of the natural universe can be seen as a three-part problem. 
Identifying the basic particles that are the constituents of matter, knowing what forces 
they feel and knowing how to calculate the behaviour of the particles given the forces. 
The Standard Model is the theory that tries to answer these questions today. 



The particles can be put into two categories, matter particles and gauge bosons. The 
matter particles, the quarks and the leptons, are the building blocks of nature. Both are 
spin- 1/2 fermions but the quarks also have colour charge as described in QCD. So far 
there are known to exist six kinds of quarks and six kinds of leptons. 

The quarks arc called up, down, strange, charm, top and bottom. They are denoted 
by their first letter and are divided into three doublets, called generations: 




The top row has electric charge q = +(2/3)e and the bottom row has q = — (l/3)e, where 
e is the magnitude of the electron's charge. 

The six leptons are also arranged in three generations: 




The electron, muon and tau have electric charge — e and each has its own neutrino of 
electric charge zero. The quarks and the leptons are the basic particles of matter. 



The second type of particles are the gauge bosons, which all have integer spin. They 
are best described together with the forces. Classically, a force is seen as a field, e.g. a 
charged object is surrounded by an electric field. But in Quantum Field Theory, which 
is based on special relativity and quantum mechanics, an important part is that particles 
and fields are treated in the same way. Whether they behave as particles or fields depend 
on the circumstances and the gauge bosons are just those particles that correspond to the 
forces and their fields. There are four forces in nature and each has its own gauge bosons: 

Force Acts on Gauge bosons 

gravitation all particles graviton 



electromagnetic all electrically 7 
charged particles 

weak quarks, leptons W"^, W~, Z° 

electroweak gauge bosons 

strong all coloured particles gj, i = 1 . . .8 

The gravitational force acts upon all particles but the strength depends on the mass and 
it can be disregarded in all but the highest energy calculations in particle physics. The 
electromagnetic force affects charged objects and holds e.g. the electrons and the nuclei 
together in atoms. It is transmitted by the photons. The weak interaction, or force, is 
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needed to explain e.g. the (3-decay of nuclei and it is transmitted by three very heavy 
bosons, which makes it a short-range interaction according to Heisenbergs uncertainty 
relation. The strong force holds the quarks together inside the nucleons as well as the 
protons and neutrons inside nuclei. The carriers are the eight gluons which themselves 
have colour charge. 

Furthermore, all particles, fermions and bosons, have antiparticles. They have 
opposite values of all charges, but the same mass, fl], || 

Particles with colour charge like the quarks and gluons can't exist on their own. They 
are permanently bound inside colourless hadrons. Hadrons can be baryons, in which case 
they are made up of three quarks, or mesons, if they are made up of one quark and one 
anti-quark. Depending on how they behave under parity-transformations and rotations, 
the mesons can be scalar-mesons, vector-mesons etc. We study properties of the three 
pseudo-scalar mesons 7r°, rj and 1]' in this article. 

1.1 Form-factors 

An interesting thing to keep in mind is that the known universe is made up of u-quarks, 
d-quarks, electrons, and the gauge bosons. All the other quarks and leptons existed at 
an early stage of the universe and can today only be seen in cosmic rays and accelerator 
experiments. Such experiments makes it possible to study the properties of the different 
interactions, i.e. how the particles behave given the forces. 

For example, studying 77 scattering provides a unique opportunity to learn the 
properties of strong interactions. Although in 77 scattering the probe and the target are 
both photons that are carriers of the electromagnetic force, they can produce a pair of 
quarks that interact strongly and are observed in the form of hadrons, e.g. pseudo-scalar 
mesons. However, the transition between a meson and two photons can't be calculated 
from QCD directly. Therefore, a temporary solution is to calculate the process without 
including all effects of the quarks and gluons and at the end take the neglected effects 
into account with an extra factor. This factor is, for historical reasons, known as a form- 
factor. The shape of this form-factor is however not known exactly. There exist different 
educated guesses which fits certain demands from QCD, but since one doesn't know 
exactly which one to use it is important to know which effects a specific choice will have. 
This article studies the variation of some physical quantities with the choice of form-factor. 

The article is organized as follows. In section 2 the form-factors are presented in more 
detail. Section 3 includes a comparison between the decay-rates for different form-factors 
for processes involving rj. Section 4 treats the cross-section of e"'"e~ e~^e~PS, where 
PS is 7r°, 7] or t]', and how it is affected by the choice of form-factor. Section 5 contains a 
discussion of the effect the form-factors have on the theoretic value of the muon anomalous 
magnetic moment and section 6 is a summary of this article. 
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2 Form-factors in Meson-Photon-Photon Transitions 



The transition between a meson and two photons is very hard to calculate in QCD. That 
depends on the fact that for low energies a^, the strong coupling constant, is too large for 
a perturbative approach to work. Beside perturbative calculations no good method exist 
in this case so the problem gets complicated and approximations have to be made. 

In the approximation that the quarks are massless, QCD has chiral symmetry. That 
means that all quarks, right-handed as well as left-handed, are treated in the same way 
and give the same result^. According to Noether's theorem, which states; 
For a system described by a Lagrangian, any continuous symmetry which leaves invariant 
the action, /£ dt, leads to the existence of a conserved current, 

this implies that the current, 5*^, which corresponds to the chiral symmetry, fulfills the 
relation 

d^S^ = 0. 

But when one combines QCD and electromagnetism, as in the vertex we are studying, 
this is no longer true. The four- divergence of the current isn't zero anymore but instead 
it has the form 

d,S^ = Ce' e^^'^^ F;, F2p, 

where C is a known constant, e the magnitude of the electron's charge, £^^"°'l^ the totally 
antisymmetric Levi-Civita tensor and F"' the electromagnetic field strength. 

This was very surprising to the people that first calculated it so it was given the name 
'the anomaly'. But since this effect exists the effective Lagrangian that should describe 
this vertex should reproduce it. That puts strong constraints on how the Lagrangian can 
look and there is actually only one way to construct it so that it fulfills this anomaly 
condition. From this Lagrangian, the term which describes the low-energy processes that 
include 7r°, rj or rj' and two photons can be written 

C = d,A,i^) A^{x) dp<j>{x) 

where is the pion decay constant, A^{x) and Aa{x) the photon fields and (j){x) the 
pseudo-scalar field. 

However, this Lagrangian only describes the process in the limit where the pseudo- 
scalar is seen as a point-like particle, with no inner structure, and the reaction is seen as 
in Fig. |1|. In reality, these particles are not pointlike but instead made up of quarks and 
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7r°, ?7, r]' 
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Figure 1: Feynman diagram for the transition when the reaction is seen as pointlike. 



^In reality, this symmetry is spontaneously broken | 
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Figure 2: The realistic situation with more quark and gluon effects. 



gluons and the reaction is more hke in Fig. 0, where the shaded area includes different 
effects from quarks and gluons. If these effects are taken into account other terms come 
in which are added to the point-like Lagrangian. The new terms all have to fulfill the 
condition d^^S'^ = 0, which means that they don't add anything to the four- divergence 
and the total Lagrangian still meets the requirement of the anomaly condition. But these 
new quark and gluon terms are very difficult to handle so they are omitted and their 
effects are only introduced as a form-factor. That is, a factor is multiplied to the vertex 
which takes into account all the effects that the extra terms would add. But since these 
effects are not completely known, the shape of the form-factor is uncertain. Still, there 
exist some constraints on it. 

2.1 Constraints on the Form-factors 

The form-factor connects three particles and therefore depends on three variables, the 
momentum-squared of the photons and the pseudo scalar. However, since we work under 
the approximation that the quark-mass is zero, the pseudo scalar will also become mass- 
less. Unlike the photons, which can be virtual, this will mean that qpg will always be zero 
and we set 

F{klklO) = Fiklkl), 
where ki and k2 are the four-momenta of the two photons. 

It would be very hard to find a form-factor that describes this transition well if one 
didn't have any constraints on how it should behave. Luckily, it is possible to calculate, 
directly from QCD, how it looks in certain limits. There are actually four different 
requirements on the shape of the form-factor. Three which are demands from QCD- 
calculations and one of a more argumentative nature. 

The first limit in which the form-factor is known is when ki = ki = 0, that is, when 
the two photons are real. Then the transition can be compared to other processes like 
7r+ fi^ u^, and it is actually in this limit that the Lagrangian is defined and is set. 
So here, by definition, all quark and gluon effects are included in the pointlike Lagrangian 
and the form-factor must fulfill the relation 

F(0,0) = L 

The next case where it is possible to calculate the result from QCD is when both photons 
are very much off shell and have the same virtuality, that is kf = k^ << 0. In this case 
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Figure 3: The high momentum imphes that the quarks move parallel to each other. 

one can use what is called the operator product expansion technique and calculate the 
matrix-element related to this transition up to order l/k"^. This gives the answer 

nklkl) = -^-^^ kl = kl = e«0, 

where Nc is the number of colours in the Standard Model, which is three. 

The third limit is when one of the photons is real and the other one very much off 
shell, e.g. kf = and << 0. Which one is zero doesn't matter since the form-factor 
has to be symmetric in the two photons. In this case one can use that this is an exclusive 
process at high momentum transfer and the transition can be seen like in Fig. ^ where 
the quark and anti-quark that builds up the meson move parallel to each other. This 
simplifies the calculations and the answer from QCD is 

F{o, e) = - ^^y , ki = 0, ki = e « 0. 

where C is a constant that is 1, 2 or 3 depending on further physical assumptions 

That F(/c^, /c|) goes like l/k"^ in these two cases is very reliable. The next terms are 
of the order l/fc^ so they are obviously much smaller. The coefficients in front are a bit 
more uncertain though, since corrections to them are of order log(fc^) and that is not so 
small. This will be referred to in the next section. 

The last constraint is not of the same form as the others. It is just that one believes 
that the quark and gluon effects which are not included in the pointlike treatment should 
somehow correspond to intermediate states with other mesons, like p. This idea comes 
among other things from the fact that the similar reaction e~^e~ 7* ^ tt+tt" is very 
well explained in a picture with an intermediate state like e~^e~ ^ 7* ^ p ^ tt+tt^. 
This means that the shape of the form-factor should be such that it can be explained 
within the picture of an intermediate state vector meson. This is not absolutely certain 
since one doesn't really know the QCD-calculations, but there are some other theoretical 
arguments that point in this direction as well [^. 

2.2 The Form-factors in this Article 

It proves very hard to find a form-factor that satisfies these four constraints. Actually 
none of the factors we study in this article does that. We have chosen four different 
form-factors, each satisfying some of the constraints, but not all. Nevertheless, it is still 
possible to see how much the choice of form-factor influences different physical quantities. 
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The first form-factor is just F{k\, = 1. Tliat means no form-factor, and the reaction 
is seen as pointlike. This form-factor obviously satisfies the first constraint but not the 
others. This is no realistic form but it can be studied to see how much the introduction 
of a form-factor means. 

The next form-factor is constructed to reproduce the experimental data that exist in 
this area. Experiments has been done with one photon real and the other off shell 0. All 
our form-factors will be compared to these data in the next section. This form-factor is 



where rUp is the mass of the p-meson. It satisfies the first and fourth requirement and with 
= 92.4 MeV and mp = 770 MeV it satisfies the third constraint within the uncertainty 
mentioned above. 

The third form- factor we have made up ourselves to satisfy three of the four constraints 
and it looks like 



It satisfies constraint one, two and three with the same comment as the last one, but not 
four, it can't be explained in terms of intermediate vector mesons. 

The last one is also constructed to satisfy a specific requirement, but this time number 
two. It is 



F{kl kl) 




{ml - kl){ml - klY 



F{kl kl) 




{ml -kl- kl)' 



F{klkl) 




{kf + k^ 




and it satisfies constraint one, two and four. 



To sum up, the form-factors we are looking at in this article are: 



Form Factor 1 : F{kl, = 1 



Form Factor 2 : F{kl,kl) 




{ml - kl){ml - kl) 



Form Factor 3 : F{k 





{ml -kl- kl) 



Form Factor 4 : F{kl,kl) 



- ^ (kf + kl 
{ml - kl){ml - kl) 
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Figure 4: Comparison between our four form-factors and data from the CLEO collabora- 
tion. 

2.3 Comparison with Experiment 

Experiments with form-factors has been done by the CLEO collaboration . They have 
studied the pion form-factor in the case that one of the photons is real and the other 
one is virtual. In this case our Form 2 and Form 3 have the same shape. In Fig. ^ the 
form-factors presented in the last section are compared with the data from the experiment. 
Their normalization of the form-factors is not the same as ours. Therefore our form-factors 
from the last section are multiplied with a factor 



647r r(7r'^ 77) 
\ (47ra)2 m^o 

It is easy to see that Form 1 and Form 4 do not satisfy the data as well as the other 
two. But still we will study the results from them to get a feeling for how much influence 
an incorrect form factor will have on various results. 
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3 Decay Rates 



The decay rate of a particle is related to its lifetime. One version of Heisenbergs uncer- 
tainty relation can be written 

AE At ~ n, (1) 

where At can be seen as the lifetime of a particle. If a particle has a finite lifetime, Eq. (|1|) 
states that it must have an uncertainty in energy. This uncertainty is usually called the 
decay width, or the decay rate, of the particle. 

In this section the decay rates of three processes, which include the meson-photon- 
photon vertex r] — >■ 77, is calculated. The reason for choosing 77 is that an upcoming 
experiment, WASA at CELSIUS, Uppsala, is going to produce up to 10^° //-particles. 
Therefore, it is interesting to conclude whether they will be able to see the difference 
between the form-factors. We will present results for our four form-factors to see whether 
they can be distinguished in the experiment. But first we study 7r° decay to check our 
approximat ions . 



3.1 77 

First we will study this process where 7r° decays. There is no dependence on form-factors 
since both photons are real, but it is still interesting to study it to check the different 
approximations we have made. The answer can be compared with the experimental 
lifetime of the pion and thus we have a measure of how well our calculations work. 



3.1.1 The Amplitude 

The process, as can be seen in Fig. ^, includes only the triangle vertex itself and no 
propagators. The triangle vertex is described by the Lagrangian 

C = d,A,i^) A^{x) dp<j>{x) e^^"'^. (2) 

Here is the proper mixture of the meson fields, 

1 12 

where e.g. the pion field is written as 

7rO(x) = / NAq) (e-'^'a, + e'^'aj) d\. (4) 




Figure 5: The decay of vr*^ into two photons. 
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The operators Ug and annihilate and create, respectively, a pion with momentum q and 
A^7ro(?) is a numerical factor which is taken care of by the Feynman rules. The other two 

meson fields can be written in the same way. 

The photon field needs an extra index since it describes a vector particle and not a 
scalar. This is introduced as a polarization vector e, which has four components. The 
photon field is then 

A^(x) = I NA(q) (e^e-'^'ak + e^'^al) d'q. (5) 

In the process we are looking at one pion is annihilated and two photons are created. The 
corresponding terms are picked out from the field operators and used in the Lagrangian 
to get the Feynman amplitude for the vertex. This gives 

M^' - 77) = ^ e'^'^'^^ qp e\MAa- (6) 

But this is not the final answer. According to the Feynman rules an extra i should be 
added and furthermore, the photons are identical particles. This means that another term 
have to be added to the amplitude, in which the substitution 71 72 is made. This gives 
the result 

A{n' - 77) = ^ e^^-^ qp + ^IMAa] ■ (7) 



3.1.2 The Matrix Element Squared 

But to calculate the decay rate we need the amplitude squared, usually called the matrix 
element squared, 

\M\^ = \A\^ = A A*. (8) 

When this is calculated, one has to sum over all polarization states. This is taken care of 
by using the rule 

^^lve^*lr = -9vT, (9) 

where the sum goes over the two polarization states of the photons and is the metric 
tensor. Rules for the Levi-Civita tensor such as 

^'Tj ^ -2{g^^gf'^ - g^Pg'^^) (10) 

also have to be used. 

The matrix element includes terms of scalar products of all the momenta in the process. 
These scalar products have to be calculated and put in. The result is 

l^lV-77) = ^-^o. (11) 
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3.1.3 The decay rate 

Now the matrix element squared should be connected to the decay rate. This is done 
using the general formula 

dV = |Mp d^r^, (12) 

where M is the mass of the decaying particle, n the number of particles in the final state 
and 

.^..^.^(P-gponp^ (13) 

In this case there are two degrees of freedom and the formula just gives 

,r = ^ (14) 

where |k| is the magnitude of either of the photons four momenta, equal for both. 

In the case we are looking at |k| = m7r/2, and since the two photons are identical and 
7r° is a pseudo-scalar, dfl = 2tt. The final result is then 

. mlo 



with a = e^/47r. 



3.1.4 The Result 

Putting in the numbers 

= 92.4 MeV a = 1/137 m^o = 134.98 MeV (15) 
gives the numerical result 

r(7r° ^ 77) = 7.734 X 10"*^ MeV (16) 

Knowing how often the pion decays in this channel, it is possible to calculate its 
lifetime, r. 

r(7c°) = = 127800 MeV'^ = 8.41 x 10"^^ s 

Experiments give 

r(7rO) = 8.4 x 10"^^ s 

so our approximations seem valid. 

This is expected to be correct for the pion but wouldn't have been as exact for rj and 
7]' since they are more complicated. First of all the s quark is much heavier than the u 
and d quark, which make our approximations more severe for rj and t]'. Furthermore their 
mass is affected by gluon effects which we haven't included. 
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Figure 6: The decay of rj into photon, electron and positron. 
3.2 77 ^ 7 e~^e~ and 77 ^ 7 fi'^i2~ 

The process with the electron and positron is described in Fig. ^ and it is the first case 
in which we can study the effects of the form factors. Here a new vertex is added which 
is described by the Lagrangian 

C = -e (eYe) A^, (17) 

where e and e are the positron and electron fields and 7 stands for the 7-matrices. There 
is also a photon propagator, which according to the Feynman rules is described by —i/kf 
in this case. These two factors get multiplied to Eq. (||) which then leads to a new 
amplitude 

A = - ^J^,p ^"'^"^ I [(e7.e)6t,fci, + (e7.e)eLA;2,]. (18) 

This is then squared, which e.g. includes working out traces of the 7-matrices. The result 
was checked with the program FORM. The matrix element squared is then used in the 
kinematics formula for three body decay, coming from Eq. ([12D , 

where ml+^- = (pe+ + Pe-)^ and m^-^ = (pe- + are the two degrees of freedom we 
have chosen. These are then integrated over, dml-^ analytically and dml+^- using the 
Gaussian quadrature program DGAUSS. In the last integration different cuts on the lower 
limit are made to see if this gives more difference between the decay rates for different 
form-factors. The decay rate is then divided by 

r*(r/ ^ 77) = 1.720 X 10"^ MeV, 

where * means that it is our theoretical value of this decay rate. We use this one because 
of the arguments before that our approximations are not quite correct in the t] and rj' 
case. In this way, the same error should exist in both numerator and denominator so the 
ratio will be a more accurate number. The result is shown in Table |l|, with numerical 
errors smaller than the last digit in the quoted values. The last line is the value without 
cut for ?7 — > 7 ii~ . We have used = 547.5 MeV . 
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cut (MeV^) 


Form 1 


Form 2,3 


Form 4 





1.620 X 


10- 


-2 


1.666 X 10- 


-2 


1.657 X 


10- 


-2 


2500 


4.614 X 


10" 


-3 


5.071 X 10- 


-3 


4.980 X 


10- 


-3 


10000 


2.580 X 


10- 


-3 


3.000 X 10- 


-3 


2.916 X 


10- 


-3 


22500 


1.507 X 


10- 


-3 


1.869 X 10- 


-3 


1.796 X 


10- 


-3 


40000 


8.601 X 


10- 


-4 


1.150 X 10- 


-3 


1.091 X 


10- 


-3 


62500 


4.614 X 


10- 


-4 


6.748 X 10- 


-4 


6.308 X 


10" 


-4 


90000 


2.237 X 


10- 


-4 


3.636 X 10- 


-4 


3.342 X 


10" 


-4 


122500 


9.258 X 


10- 


-5 


1.705 X 10- 


-4 


1.539 X 


10- 


-4 


160000 


2.958 X 


10- 


-5 


6.322 X 10- 


-5 


5.585 X 


10- 


-5 


202500 


5.825 X 


10- 


-6 


1.488 X 10- 


-5 


1.284 X 


10- 


-5 


250000 


3.394 X 


10- 


-7 


1.077 X 10- 


-6 


9.050 X 


10- 


-7 




5.506 X 10" 


-4 


7.744 X 10- 


-4 


7.284 X 10" 


-4 



Table 1: Decay rate for rj ^ j e^e , normalized to r(?7 77). 



3.2.1 Experimental Uncertainties 

It is necessary here to get a feeling for how large the uncertainty will be in the forthcoming 
experiment. A first estimate of the error is V^, where N is the number of decaying 
particles. is given by 

N = No^^ = No^ BR{r^ 77), (20) 

where Fq is the total decay rate and F/F^^ the number given in the table. From this one 
can calculate the uncertainty on the normalized decay rate, 



F/F^^ 



(21) 



\ BR{7j ^ 77) No 
Without any cut and with A^o = 10^° and 

BR{r] 77) = 0.39, 

this gives an uncertainty a ~ 10^^. So it should be possible to discern the different 
form-factors here. 

But then a cut is introduced to see whether the difference between the form-factors 
will increase. For the cut on 160000 MeV"^, the difference between the form-factors is 
quite sizeable. From Eq. (|2T|) one sees that a ~ 10"'' in this case, which means that it 
should still be possible to distinguish between the three different form-factors with a high 
cut in the experiment. The reaction — 7 fi~^fi~ gives the same result. 

It remains to check whether it is possible to see the difference between Form 2 and 3 
in some other decay. 

3.3 ij fi^jj,^ e+e^ 

As one can see in Fig. |^, one new vertex and one new propagator is added to the last 
reaction. This leads to an amphtude 

• 4 

A = - ^^l^ p ^"^"^^ [inae)h,m^,fi) + (e7,e)fc2M(/i7a/i)]- (22) 
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Figure 7: The decay of 7] into muon, anti-muon, electron and positron. 

Calculating the corresponding matrix element squared includes working out the Dirac 
traces and once again the result is checked with FORM. The matrix element is then 
inserted in the formula for the four body decay rate, which has five degrees of freedom, 

= '^'ItlfJ''^'**' d{cose^*) t/(cos0/*) rfC dml dml (23) 



where mi = \J {pe+ +Ve-Yi ^2 = y {Pfi+ + Pfi-)^ *{**) indicates that the quantity 
should be calculated in the rest frame of 71(72)- These five degrees of freedom are 
again integrated over, this time the angle integrations are done analytically and the 
other two with DGAUSS. This leads to different results for the form-factors according 
to Table 0. Once again cuts are made, this time we impose a lower limit on ml, to see 
if this enhances the difference between the factors. For the lowest cut, 625 MeV"^, the 



cut (MeV^) 


Form 1 


Form 2 


Form 3 


Form 4 





3.99 X 10" 


-6 


5.61 X 10' 


-6 


5.62 X 10' 


-6 


5.28 X 10' 


-6 


625 


1.39 X 10" 


-6 


1.96 X 10' 


-6 


1.97 X 10' 


-6 


1.84 X 10' 


-6 


2500 


8.26 X 10' 


-7 


1.16 X 10' 


-6 


1.17 X 10' 


-6 


1.09 X 10' 


-6 


5625 


5.21 X 10' 


-7 


7.37 X 10' 


-7 


7.44 X 10' 


-7 


6.92 X 10' 


-7 


10000 


3.29 X 10' 


-7 


4.68 X 10' 


-7 


4.74 X 10' 


-7 


4.39 X 10 


-7 


15625 


2.03 X 10 


-7 


2.91 X 10' 


-7 


2.96 X 10 


-7 


2.73 X 10' 


-7 


22500 


1.20 X 10' 


~7 


1.74 X 10' 


-7 


1.77 X 10' 


-7 


1.63 X 10 


-7 


30625 


6.64 X 10' 


-8 


9.79 X 10' 


-8 


9.99 X 10 


-8 


9.12 X 10' 


-8 


40000 


3.36 X 10' 


~8 


5.04 X 10' 


-8 


5.16 X 10' 


-8 


4.68 X 10' 


-8 


50625 


1.49 X 10' 


-8 


2.28 X 10' 


-8 


2.34 X 10' 


-8 


2.11 X 10' 


-8 


62500 


5.36 X 10 


-9 


8.46 X 10' 


-9 


8.70 X 10' 


-9 


7.78 X 10' 


-9 


75625 


1.37 X 10' 


-9 


2.22 X 10' 


-9 


2.29 X 10' 


~9 


2.03 X 10' 


~9 


90000 


1.68 X 10" 


10 


2.82 X 10" 


10 


2.92 X 10" 


10 


2.57 X 10^ 


10 



Table 2: Decay rate for 77 — fi^fi e^e , normalized to r{r] — ^ 77). 



uncertainty calculated from Eq. (|21|) is a ~ 10 ^. Then it is possible to distinguish 
Form 1 and Form 4 from the others. For a higher cut, like 30625 MeV'^, the error should 
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be 0" ~ 5 X 10 ^. It still is not possible to see the difference between Form 2 and Form 3. 



3.4 T] e~^e e~^e 

The last decay we will study is the one depicted in Fig. ||. Kinematically, it is the same 
as the last one but the amplitude gets more complicated. A start is to take the amplitude 
from the last case, Eq. (^2j) and change /i to 62- This gives 

• 4 

^ = - wfV^ ■M:2 [(ei7aei)A;i^(e27.e2) + (ei7.ei)A;2;.(e27ae2)]. (24) 

oV ovr i^TT "'i"'2 

But this isn't enough, because here we once again have identical particles. Two pairs 
of two identical particles give four diagrams, which can be divided into pairs which give 
equal results. The diagram in Fig. || represents one of the pairs and Fig. | represents 
the other pair. Both of these diagrams have to be taken into account and according to 
the Pauli principle they come in with different signs. Furthermore, this process is second 
order and this implies, according to the Feynman rules, an extra factor of 1/2. Finally, 
the two pairs don't have the same propagators. In the first case they are 

h = Pe+ + Per 
h = Ve+ + Ve- 

but in the second case the propagators are 

h = Pe+ + Pe- 
k4=Pe+ +Pe-- 

This gives an amplitude 



with 




Figure 8: The decay of rj into two electrons and two positrons. 
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73(^3, es) 



ri{q) 




74(^4, 64 



Figure 9: The other electron configuration. 



and 



A, 



(617062) fc3M(e27t^ei) + (6171^62) fc4At (627^61^ 

3 4 



(27) 



This amphtude is then squared, using FORM, and put into Eq. ([23|), where an extra 
factor 1/4 is added since there are identical particles in the final state, to get the decay 
rate. The five degrees of freedom are integrated over using the Monte Carlo integration 
program VEGAS. Cuts are made on the low limit of kf, and k^, and the result is 

presented in Table ^. The errors in parantheses are the errors quoted from VEGAS. The 
errors on the numbers without parantheses is maximum one in the last digit. 



cut (MeV^) 


Form 1 


Form 2 


Form 3 


Form 4 





6.40(2) X 10-^ 


6.73(2) X 10-5 


6.71(2) X 10-5 


6.64(2) X 10-5 


625 


9.93 X 10-*^ 


1.10 X 10-5 


1.10 X 10-5 


1.08 X 10-5 


2500 


3.83 X 10-6 


4.40 X 10-6 


4.41 X 10-6 


4.29 X 10-6 


5625 


1.57 X 10-^ 


1.86 X 10-6 


1.86 X 10-6 


1.80 X 10-6 


10000 


6.07 X 10-^ 


7.42 X 10-^ 


7.46 X 10-^ 


7.15 X 10-^ 


15625 


2.07 X 10-^ 


2.63 X 10-^ 


2.64 X 10-^ 


2.51 X 10-^ 


22500 


5.80 X 10-* 


7.65 X 10-* 


7.73 X 10"* 


7.28 X 10-* 


30625 


1.24 X 10-* 


1.73 X 10-* 


1.76 X 10"* 


1.63 X 10-* 


40000 


1.86 X 10-9 


2.74 X 10-9 


2.79 X 10-9 


2.56 X 10-9 


50625 


1.40 X 10-1° 


2.19 X 10-1° 


2.25 X 10-1° 


2.01 X 10-1° 


62500 


1.74 X 10-12 


2.84 X 10-12 


2.96 X 10-12 


2.61 X 10-12 



Table 3: Decay rate for t] e^e e^e , normalized to r(?7 — > 77). 



The difference between Form 2 and 3 is very small until the cut is 10000 MeV"^. 
According to Eq. (|21|), the uncertainty here is cr ~ 10"*, so it is not possible to see 
the difference here. For the highest cut the distinction is larger but the experimental 
uncertainty is up to the size of the value itself. So, unfortunately, it seems not to be 
possible to choose between form-factor 2 and 3 from decay experiments with this kind 
of precision. In the next section we will see if it is possible in collision experiments like 
e~^e~ —>■ e~^e~PS, where PS is 7r°, rj or 77'. 
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4 Cross Sections 



The cross section for a reaction is related to the probabihty for it to happen. The name 
comes from the classical picture where the total cross section for e.g. pp scattering is 
the area in which both particles have to be when they pass each other for a reaction to 
happen. The cross section for e.g. pp tt^tt^ is then this area times the probability that 
precisely tt+tt^ is in the final state and nothing else. 

The process e+e" — >■ PS e+e", where PS is 7r°, rj or t]', includes a meson-photon- 
photon vertex. Therefore this is another opportunity to study the effects different form- 
factors will have. We will calculate the cross section for all three processes with our 
four different form-factors to see whether the difference will be larger here, possibly large 
enough to distinguish between Form 2 and 3. 

Experimentally, these cross sections can be studied both at CESR, US, and at 
DA$NE, Italy. At Cornell they are studied at CLEO with a center of mass energy 
EcM = 10.58 GeV and in Italy at KLOE with Ecm = 1-0192 GeV. We will look at both 
these energies in our calculations. 



4.1 The Matrix Element 

Irrespective of the pseudo-scalar created, two diagrams contribute to the process, depicted 
in Fig. and Fig. |TT]. As seen, the process includes the same vertices and propagators 
as the earlier treated process rj e~^e~ e~^e~. Then, for reasons of crossing symmetry, 
the matrix element should also look the same, with a different numerical factor for 7r° and 
Tj'. But not exactly the same. The creation of one of the electron-positron pairs before 
corresponds to annihilation of a positron-electron pair now. This taken into account, the 
amplitude is 



PS e+e- 



le 



{2i 



with 



A,, 



(ei7t.ei)fci^(e27ae2) + (ei7aei)fc2^(e27^e2) 

n'in'2 

(ei7ae2)fc3M(e27i^ei) + (ei7ve2)fc4M(e27aei) 

^-2/^2 




PS{q) 



Figure 10: Diagram contributing to the cross section for e+e PS e+e . 
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Figure 11: Second diagram contributing to the cross section for e~^e PS e^e 
and 

f 8 for 7r° 

C = < 8^3 for ri 
[ 4^372 for ri 

where the difference from the earher case hes in which part of the electron and positron 
fields that are used. In the calculations we have used an effective C for 77 and 77', 
which reproduces the correct decay rates for them, Yiji 77) = 4.6 x 10~^ MeV and 
Y(n' ^ 77) = 4.27 X 10-3 MeV. 

The amplitude is then squared and the result is the matrix element squared. 

4.2 The Kinematics 

The general formula for the cross section with three particles in the final state is 

4\/(pi ■ - mlml 



with 

,^ .4/ , N d^P^ d^P^ d^P 

d$3 = 5 (Pl + P2 - P3 - P4 - P5 



(27r)3 2Es (27r)3 2E4 (27r)3 2E5 

where pi and p2 are the four momenta of the incoming particles and pa, p4 and ps are 
the momenta of the outgoing particles. In this case (i$3 can be rewritten as 

\K+\d^l+ \Pe-\d^e- dml + 

= — T^^TTT^ (30) 

16(27r)y mpg^+ Ecu 



where rnpg^+ = ^{q + P^+Y and * means that the quantity is calculated in the rest frame 

of €2 and the meson. Furthermore, since the lab system, where we want to calculate our 
quantities, is the same as the Ecu system we have 

Pl = |^,0,0,W^-mn (31) 
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P2 = 1^^.0,0, -\/^ w? 

which imphes 



e I 5 



[Pi • V2Y - mlml = EcM \Pin\, (32) 
where is either pi or p2. 

4.3 The Result 

The z-axis is set by the incoming particles, but one can rotate freely around it. This 
means that one of the angles in dfl^- becomes just 27i. Furthermore, a factor 1/4 is added 
since one has to take the average over incoming spins. Then there are four degrees of 
freedom left and the expression for the cross section is 

bl+l \Pe-\ |Mp 

da = / -2 p-r dm + dVt*+ cos 6*2. (33) 

The matrix element squared is made up of scalar products of the four vectors 
p- and p- . Therefore one has to figure out the how these vectors look in some common 
frame. This is done by expressing them in a frame where they are known and then rotating 
and boosting them to the lab frame. 

Then everything is known and one can integrate to get the full cross section. This 
is again done by using VEGAS integration on all four degrees of freedom. Cuts are 
introduced, this time directly on the angle between the outgoing electron and positron 
and the z-axis. The lowest cut is constructed so that the electron and positron precisely 
go through the tracking system of the detectors. 

The results, in pbarns, are presented in Table ^ to Table |^. The numerical errors are 
maximum one in the last digit. 

We have used the following masses 

m^o = 134.98 MeV = 547.5 MeV m^' = 957.8 MeV, 



cut ( cos6' <) 


Form 1 


Form 2 


Form 3 


Form 4 


0.9 


2.41 


2.09 


1.94 


1.63 


0.8 


1.38 


1.33 


1.21 


0.919 


0.6 


0.596 


0.690 


0.627 


0.409 


0.4 


0.266 


0.355 


0.328 


0.191 



Table 4: Cross section for e~^e ir e~^e with Ecm = 1.019 GeV and extra cut. 
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cut (1 cos^ 


<) 


Form 1 


Form 2 


Form 3 


Form 4 


0.93 
0.8 
0.6 
0.4 


4.43 
1.76 
0.762 
0.355 


6.83 X 10"^ 
3.69 X 10"^ 
1.23 X 10-5 
6.85 X 10-6 


6.20 X 10-=^ 
5.25 X 10-^ 
8.58 X 10-5 
2.39 X 10-5 


3.56 X 10-2 
2.05 X 10-2 
1.32 X 10-2 
8.04 X 10-3 



Table 5: Cross section for e+e 7r° e+e with Ecm — 10.58 GeV and extra cut. 



cut ( 


cos^ 


<) 


Form 1 


Form 2 


Form 3 


Form 4 


0.9 
0.8 
0.6 
0.4 


1.26 
0.694 
0.291 
0.123 


0.974 
0.557 
0.260 
0.126 


0.963 
0.548 
0.255 
0.123 


0.871 
0.457 
0.186 
0.0800 



Table 6: Cross section for e+e — > rj e+e with Ecm ~ 1.019 GeV 



cut ( 


cos 6 


<) 


Form 1 


Form 2 


Form 3 


Form 4 


0.93 


3.90 


6.06 X 10-^ 


5.48 X 10-3 


3.15 X 10-2 


0.8 


1.55 


3.27 X 10-5 


4.65 X 10-^ 


1.81 X 10-2 


0.6 


0.673 


1.09 X 10-5 


7.63 X 10-5 


1.17 X 10-2 


0.4 


0.314 


6.05 X 10-6 


2.12 X 10-5 


7.10 X 10-3 



Table 7: Cross section for e+e rj e+e with Ecm — 10.58 GeV and extra cut. 



4.3.1 Additional Cuts 

As stated earlier, some of the form-factors corresponds to a physical picture with interme- 
diate mesons e.g. p. Therefore the mathematical expression of those factors are related 
to a meson propagator, 

where k is the four momenta and M the mass of the meson. This is however only an ap- 
proximation of the full propagator, which should include different higher order corrections, 
such as loops. This would have changed the propagator to 

(35) 



A;2-M2-E' 



where E is a complex quantity related to the decay rate of the propagating meson. 

Using this approximation introduces problems when — M'^, which will happen in 
some cases with the energies above. In the case with the high energy, Ecm = 10.58 GeV, 
the problem arises with all three of the particles. This is solved by introducing a new cut. 
AU events in which kj, kj, kj or kj is in the range {Mp - 300 MeVf (Mp + 300 MeVf 
are not counted in the cross section. The reason for choosing 300 MeV is that it is twice 
the decay width of the p. A good approximation for E is iMpFp, which means that when 
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cut ( 


cos^ 


<) 


Form 1 


Form 2 


Form 3 


Form 4 


0.9 
0.8 
0.6 
0.4 


0.205 
0.114 
0.0457 
0.0172 


0.183 
0.100 
0.0391 
0.0145 


0.183 
0.100 
0.0392 
0.0146 


0.187 
0.102 
0.0399 
0.0148 



Table 8: Cross section for e+e rj' e+e with Ecm — 1.019 GeV 



cut ( cos^ <) 


Form 1 


Form 2 


Form 3 


Form 4 


0.93 


6.61 


1.05 X lO--'^ 


9.38 X 10-3 


5.38 X 10-2 


0.8 


2.65 


5.65 X 10-5 


8.02 X 10-^ 


3.09 X 10-2 


0.6 


1.15 


1.87 X 10-5 


1.32 X 10-4 


2.00 X 10-2 


0.1 


0.536 


1.03 X IQ-'' 


3.65 X 10-' 


1.21 X 10-2 



Table 9: Cross section for e'^e — > rj' e^e with Ecm = 10.58 GeV and extra cut. 



q — Mp±2rp, the absolute value of the term we have neglected can at most be one fourth 
of q'^, and therefore should not make a large difference. 

The same thing happens with the low energy and the pion. But this time it is not 
possible to use the same cut as last time, since all our events lie in that region. We solve 
this by cutting in a more narrow interval, Mp ± 200 MeV. Then the neglected term can 
be up to half the size of the terms left but that will have to do. 

Experimentally, it is all right to make a more narrow cut in the low energy case. A 
narrow cut means that a small error in the measurement of the energy will make much 
difference, but this is compensated by the fact that low energy means better resolution, 
experimentally. 

4.3.2 Experimental Uncertainties 

Once again it is interesting to calculate the uncertainties in the experiments. The number 
of events is in this case given by the formula 

N = La At, (36) 

where L is the luminosity of the experiment and At the running time. In both experiments 
the standard luminosity should be L = 5 x 10^2 cm~^ and we will calculate with one 
year of running time. At = 10^ s. 

As mentioned before, a first estimate of the uncertainty is ±\/]V, where is the 
number of events. From this one can get 

which is the result we use. 
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4.3.3 Conclusions 



Just as in the decay rates it is possible to distinguish between Form 1 and 4 in all of these 
experiments. However, here the difference is large between Form 2 and 3 in some cases 
as well. 

In the high energy experiment it should actually be possible to discern all four form- 
factors for all three particles. For e.g. rj with the cut |cos^| < 0.93, the results with 
errors are 

Form 1 (3.90 ± 0.03) 

Form 2 (6 ± 3) xlO"^ . . 

Form 3 (5 ± 1) xlO'^ ^ ^ 

Form 4 (3.2 ± 0.3) xlO'^ 

In the low energy experiments it is more difficult. There it is only in the case with the 
pion that it should be possible to distinguish between all four. The results, with cut 
I cos^l < 0.9 are 

Form 1 (2.41 ± 0.02) 

Form 2 (2.09 ± 0.02) . . 

Form 3 (1.94 ±0.02) ^ ^ 

Form 4 (1.63 ±0.02) 

On the other hand, the difference is visible for all the four different cuts in this case. 

The above numbers are calculated with the luminosity mentioned above and a typical 
year, which totals to an integrated luminosity of 5000 pb~^. 
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5 The Muon Anomalous Magnetic Moment 



A magnetic moment is associated with a rotating charged body. Classically this magnetic 
moment fi is related to the angular momentum I of the rotation through 

(40) 

where q is the charge and m is the mass of the body. But particles also have intrinsic 
angular momentum, spin. This also gives rise to a magnetic moment according to 

f^s = gs^ s, (41) 

where s is the spin angular momentum vector and q and m is the charge and the mass of 
the particle. Qs, known as the gyromagnetic ratio, is just a constant which, according to 
Dirac theory, is exactly equal to 2 for electrons and muons. 

However, in Quantum Field Theory, which treats particles as quantized fields as men- 
tioned in the introduction, this is no longer true. Fig. ^ depicts the lowest order diagram 
for the interaction between a muon and a magnetic field, which gives gs = 2. 




Figure 12: Lowest order diagram for muon scattering by a magnetic field. 

But then there are higher order corrections, which will change the amplitude for the 
interaction between the muon and the magnetic field. This can be taken care of by chang- 
ing the value of the gyromagnetic ratio by a small amount. To study these corrections 
one defines the anomalous magnetic moment, a^, as 

a, = ^igs-2). (42) 

This will then be a measure of how large the higher order corrections from Quantum Field 
Theory are. At present the theoretical standard model estimate for is 0, H, ^ 

a^{th) = 11659178(7) x 10~^°, 

where the estimated error is in parentheses, gets contributions from five different types 
of diagrams. 

1. Pure QED contributions are the largest corrections. They consist e.g. of extra 
photon propagators or quark- and lepton loops, see Fig. |l^. 
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Figure 13: Higher order QED diagrams for muon scattering by a magnetic field. 



2. Hadronic vacuum polarization, which means that all kinds of hadron effects are 
added to the photons, e.g. quark-loops or intermediate mesons. 

3. Higher order hadronic vacuum polarization effects, which means even more hadron 
effects. 

4. Hadronic light-by- light scattering contributions. These are more complicated split- 
tings and loops in the photon from the magnetic field. The dominant diagram is 
depicted in Fig. 0. 

5. Electroweak contributions, which means that Z^- and TV^-propagators are added 
in different ways. 



The piece that is studied in this paper is the one depicted in Fig. pi since it includes two 



meson-photon-photon vertices and therefore also depends on the form-factors. We will 
calculate the contribution from this process to the anomalous magnetic moment for our 
form-factors 2,3 and 4 and compare the results. 

A running experiment at Brookhaven National Laboratory measures the anomalous 
magnetic moment for the muon. In the future they hope to get an accuracy around 
±4 X 10~^°, improving by a factor more than twenty the previous experimental determi- 
nation at CERN. Therefore, it is important that the theoretical uncertainties are lowered 




Figure 14: The dominant light-by-light contribution to the muon magnetic moment. 
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to the same order as the aimed BNL uncertainty. Then the measurement can become 
a precision test of the quantum corrections of the electroweak sector of the Standard 
Model. Even more, with the same low theoretical and experimental uncertainties, and 
when combined with other high precision results, could become an excellent probe of 
physics beyond the Standard Model. 



5.1 The Result 

All of the hadronic light-by- light contributions can be depicted as in Fig. where the 
shaded blob includes all kinds of hadronic effects. The numbering in the figure corresponds 
to one of six permutations, due to the photons being identical particles. The amplitude 
for this process is 

A = I |7 4 / d% f d% 1 

^ J (27r)4 J (27r)4 q^plp^pl - m^){pl - m^) ^ ' 

y,W°'l^{p^^p2,p^)u{p')'^a{h + ^)lu{h + m)-fp u{p) 



-|-five more permutations, 

where YiP'^°'l^ ^ called the four-point function, describes the shaded blob. In our case, Yi''^'^^ 
is just the product of the vertices and propagators corresponding to Fig. For low 
momentum transfer, this amplitude can be expanded in a Taylor series in {p — p'). The 
linear term is then, with the momentum transfer set to zero, by definition the anomalous 
magnetic moment. We calculate the pseudo-scalar light-by-light contribution to this term 
for form-factor 2, 3 and 4. In doing that, the integrations in Eq. ( ^5] ) should go from 
zero to infinity. But it is interesting to see from which energy-regions the contributions 
come. Therefore, we introduce a cut on the upper limit. This cut is then raised until the 
result converges to some number, which means that a^^ has been saturated. The result is 
presented in Table |10| when PS = vr". Table [Tl| when PS = rj and Table when PS = rj'. 
The numbers in parentheses are the errors quoted by VEGAS. The contributions from 
the three pseudo scalar mesons should then be added to get the final value of a^^. The 




p' P4 P5 P 



Figure 15: Hadronic light-by-light contribution to a^. 
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cut 


Form 2 


Form 3 


Form 4 


(GeV) 


(xlO-iO) 


(xlO-iO) 


(xlO-iO) 


0.4 


-2.683(2) 


-2.707(2) 


-2.821(2) 


0.5 


-3.440(3) 


-3.449(3) 


-3.683(3) 


0.7 


-4.448(5) 


-4.622(5) 


-4.936(5) 


1.0 


-5.152(7) 


-5.538(8) 


-5.972(8) 


2.0 


-5.60(1) 


-6.43(1) 


-6.96(1) 


4.0 


-5.64(1) 


-6.67(2) 


-7.21(2) 


8.0 


-5.63(1) 


-6.74(2) 


-7.27(2) 


12.0 


-5.63(1) 


-6.70(2) 


-7.25(2) 



Table 10: a^^ when PS = 7r°. 



cut 


Form 2 


Form 3 


Form 4 


(GeV) 


(xlO-i°) 


(xlO~i°) 


(xlO-i°) 


0.4 


-0.4029(2) 


-0.4071(2) 


-0.4289(2) 


0.5 


-0.5819(3) 


-0.5936(3) 


-0.6353(3) 


0.7 


-0.8678(6) 


-0.9088(6) 


-0.9967(6) 


1.0 


-1.1058(9) 


-1.211(1) 


-1.357(1) 


2.0 


-1.281(1) 


-1.545(2) 


-1.768(2) 


4.0 


-1.301(2) 


-1.647(3) 


-1.900(3) 


8.0 


-1.302(2) 


-1.670(3) 


-1.937(3) 


12.0 


-1.301(2) 


-1.676(4) 


-1.943(3) 



Table 11: when PS = rj. 



rcsu-ltj IS 

Form Factor 2 -7.94 x 10~^° 

Form Factor 3 -9.73 x 10"^° (44) 

Form Factor 4 -10.85 x 10"^° 

If nothing else were known of the form-factors, one would have to take the mean value of 
all these and an uncertainty that covered all of them. Then the final result would be 



= -(9.51 ± 1.6) X 10"^°. (45) 



However, it looks like it is possible to exclude Form 4 already from the experiment men- 
tioned in Section 2.3. Therefore, we take a mean value of Form 2 and 3, which gives 

aj"^ = -(8.83±0.9) X 10"^°. (46) 

This is our final result for a^^ and since it isn't possible to discern Form 2 and 3 yet, this 
must be the smallest possible uncertainty on this value. This can be compared with the 
result in reference P 

= -(8.5 ± 1.3) X 10-1° (47) 

and m 



a^s = -(8.27 ± 0.64) x 10"^° (45 
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cut 


Form 2 


Form 3 


Form 4 


(GeV) 


(xlO-iO) 


(x 10-10) 


(xlO-iO) 


0.4 


-0.2595(1) 


-0.2623(1) 


-0.2770(1) 


0.5 


-0.3882(2) 


-0.3964(2) 


-0.4260(2) 


0.7 


-0.6113(4) 


-0.6417(4) 


-0.7102(4) 


1.0 


-0.8168(6) 


-0.9008(6) 


-1.0262(7) 


2.0 


-0.9874(9) 


-1.219(1) 


-1.443(1) 


4.0 


-1.008(1) 


-1.324(2) 


-1.597(2) 


8.0 


-1.010(1) 


-1.351(2) 


-1.645(2) 


12.0 


-1.009(1) 


-1.353(2) 


-1.654(3) 



Table 12: a^^ when PS = 77'. 



Even if neither Form 2 nor Form 3 satisfies all the physical demands on a form-factor, 
our result means that one can not use arguments from existing experiments to lower the 
uncertainty more than ours. To get better precision, experiments like the cross sections 
mentioned in Section 4 have to be done, which can better discern between different form- 
factors. 

6 Summary 

In this thesis we have been studying form-factors in meson-photon-photon transitions. 
In this kind of transitions the exact form of the form factors is not known and 
therefore different types are investigated. We have studied four different factors, which 
all satisfy some, but not all, of the constraints on form-factors coming from QCD. They are 

Form Factor 1: F(A;?,/c|) 
Form Factor 2: F{kl,kl) 

Form Factor 3 : F{kj, k^) 



Form Factor 4 : F{kj, /c^) 



The idea is to calculate a set of physical quantities whose value depend on the choice 
of form-factor. This is done to see if it should be possible to distinguish between the 
four in experiment. Already in section 2, results from an experiment done by the CLEO 
collaboration were presented, which more or less rules out Form 1 and 4. 

In section 3 the decay rates of processes involving — 77 are calculated. The results 
are compared with the expected uncertainty from an upcoming experiment WASA at 
CELSIUS, Uppsala. Here again it seems like they will be able to exclude Form 1 and 4 
but not choose between Form 2 and 3. 



= 1 



{ml-kl){ml-k^] 
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In section 4 the cross section of three processes, which include a meson-photon-photon 
transition, are calculated. The results are again compared with estimated uncertainties, 
this time from experiments at CESR, US and at DA$NE, Italy. It turns out that in 
this type of experiments it should be possible to make a more precise measurement of 
form-factors and possibly discern all four. 

Finally, in section 5, we make some comments about how the present uncertainty 
around the choice of form-factor affects the value for the muon anomalous magnetic mo- 
ment. When examining the pseudo-scalar part of the light- by- light contribution to a^, 
we can estimate the lowest possible error in the theoretical value by looking at how much 
it differs when using Form 2 and 3. Since no present experiment can discern between 
the two, the uncertainty must be at least large enough to cover both values. To get a 
better precision, experiments have to be done, which can better discern between different 
form-factors. This is important since the muon anomalous magnetic moment is one of the 
most precisely measured quantities in physics and furthermore can be a probe for physics 
beyond the Standard Model. 
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